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Abstract 

A calculation is presented of the radiative decay of the T(nS) into a bound state of bottom 
squarks. Predictions are provided of the branching fraction as a function of the masses of the 
bottom squark and the gluino. Branching fractions as large as several times 10~^ are obtained for 
supersymmetric particle masses in the range suggested by the analysis of bottom-quark production 
cross sections. Data are shown that limit the range of allowed masses. Forthcoming high-statistics 
data from the CLEO Collaboration offer possibilities of discovery or significant new bounds on the 
existence and masses of supersymmetric particles. 
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A long-standing puzzle in high-energy strong interactions is the fact that the rate for bot- 
tom quark {bb) production at the Fermilab Tevatron collider is two to three times greater 
than the theoretical prediction from quantum chromo dynamics (QCD) The next-to- 
leading order QCD contributions are large, and a combination of further higher-order effects 
in production and/or fragmentation may eventually reduce the discrepancy [2]. An alterna- 
tive explanation is offered in Ref. j^l in the context of physics beyond the standard model. 
There, it is argued that a solution may be provided by the existence of a light bottom 
squark b and a light gluino g, with masses in the ranges 2 GeV < < 5.5 GeV and 
12 GeV < nig < 16 GeV. The g and the b are the spin-1/2 and spin-0 supersymmetric 
partners of the gluon (g) and bottom quark (6). The masses of all other supersymmetric 
particles are assumed to be arbitrarily heavy, i.e., of order the electroweak scale or greater P|. 
While this scenario is not among the more popular schemes for supersymmetry breaking, 
the hypothesis of a light bottom squark is not inconsistent with direct experimental searches 
and indirect constraints from other observables 0,0, 0,0,0] • Therefore, it is essential either 
to confirm or to refute this proposal by examining its implications for additional processes. 
In this Letter, we show that high-statistics data being accumulated and analyzed now at 
the Cornell Electron Storage Ring (CESR) facility [9J could provide definitive confirmation 
of the proposal of Ref. |3| or severely constrain the allowed parameter space . 

In the proposal of Ref. 0], it is possible that the bottom squark is relatively stable and, 
hence, bound states of a bottom squark and bottom antisquark (sbottomonium) could exist. 
These bound states could be produced in radiative decays of bottomonium states, such as 
T 5*7, where 5* is the S-wave bound state of a bb* pair. Alternatively, the b could decay 
promptly via i?-parity and baryon- number violation P, Q], and no bound state would be 
formed. 

In this Letter, we compute the rate for an T{nS) to decay radiatively into an S'-wave bb* 
bound state. We show that, provided that a bound state is formed, the resonance search 



by the CUSB Collaboration 



ll| already increases the allowed lower bounds on and rrig. 



Discovery of the bound states may be possible with the high-statistics 2002 CLEO-c data 
set, or a larger range of bottom-squark and gluino masses may be disfavored. 

The mass eigenstates of the bottom squarks, bi and 62, are mixtures of bi and the 
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FIG. 1: Feynman diagrams for the process bb bb*'j. 



supersymmetry partners the left-handed (L) and right-handed (R) bottom quarks: 

|5i) = sin6i\bL) +cos6i\bR), (la) 
\b2) = cos 9~,\bL) -sin e-,\bR). (lb) 

We take bi to be the eigenstate of hghter mass, and we drop the subscript 1 in the remainder 
of this paper. The mixing angle is constrained by the requne-nent that the couphng of a 

bb* pair to the Z boson be sufficiently small to be compatible with data At lowest order 
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(tree-level), this requirement implies that sin^ ^^j; ^1/6 

We calculate the decay rate r(T — >■ 5*7) in the framework of the nonrelativistic QCD 
(NRQCD) factorization formalism ^ First, we compute the amplitude in full QCD 

for the process b{p) + b{p) —>■ b{q) + b*{q) + '~f{k). Typical Feynman diagrams are shown 
in Fig. m The indices i, j, k, and / label the colors of the incident b and b and the final b 
and b*, respectively. The color index of the exchanged gluino is a. The Feynman rules [lol | 
depend on the mixing angle O^. 

We carry out the computation in the bb rest frame and choose the radiation gauge. Then, 
the photon-6-6* vertex vanishes, since it is proportional to e* ■ (2g + k) = e* ■ (2p), where e 
is the polarization of the photon. Therefore, we need to compute only the two diagrams of 
the type in Fig. ^a), in which the photon attaches to the b or the b. We make use of the 



^ In the first paper of Ref. 12|, it is argued that one- loop contributions may render the hght g and hght b 
scenario inconsistent with data at the 2a level, unless the mass of the heavier 62 is less than about 125 
GeV. Making somwhat different assumptions, the authors of the second paper obtain a 5(j bound of 180 
GeV and the author of the third paper obtains a 3cr bound of more than 200 GeV. The possibility that 
the mass of the &2 could be as low as 100 GeV or so is not excluded by data. Analysis of data from e~^e~ 
interactions at LEP, for example, would require first a detailed modeling of the decay modes of the 62. 

^ One can adapt the NRQCD formalism for spin-1/2 quarks to the case of spin-0 squarks by dropping the 
spin-dependent interactions and replacing Pauli fields with scalar fields. 



spin-triplet projector [l^ 

J2 «(P> ^iM-p, A2)(i Ai i A2II er) = --^ /r{^ - m,), (2) 

Ai,A2 

where ex is the polarization of the bb system, the Dirac spinors u{p,si) and v{p,S2) 
have the relativistic normalization u{p,s)u{p,s) = v{p, s)v{p, s) = 2Ei,{p), and Ei{p) = 



y'pf+rnf. Projecting onto color-singlet bb and bb* states, which leads to the color factor 
{6ij6ki/N^)T^^T!li = |, we find that the amplitude in full QCD is 

M. = i — ^^^giech sin Qi cos d-, e ■ ex- (3) 
6mg 

We use plane- wave states for the quarks, squarks, and photon, without any factors \/{2E). 
Squaring the matrix element and averaging over the polarizations of the T, we obtain 

\W=^9yel^^'e,cos-e,. (4) 

We match the squared matrix element in Eq. onto NRQCD for the bb system in the 
hb rest frame and onto NRQCD for the bb* system in the bb* rest frame using, respectively, 



= (5) 



\M? = F,CSo)mOfCSom. (6) 

Here, Fi and Fi are short- distance coefficients, Oi{^Si) is a four-quark operator, and 
^i''*(^'S'o) is a four-squark operator: 

0,CS,) = ^^a'xX^a'^P, (7) 
ofCSo) = ri'J2\~^b*+X){bb* + X\rx. (8) 

X 

The Pauli field ip annihilates a b quark, the Pauli field x creates a b antiquark, the scalar 
field ip annihilates a b squark, the scalar field x creates a b squark, and the sum is over 
all intermediate states X. The short-distance coefficient Fi{^Si) for the bb operator con- 
tains {0\Of'* Sq)\0) , and the short-distance coefficient for the bb* operator Fi(^S'o) contains 
\Oi{^Si)\bb) . In the Born approximation, the matrix elements are 



OiCSi)\bb) = 2{2EhfN,, (9) 
{0\dfCSom = {2E-,fN,. (10) 



where Ej^ is the energy of b{b*) in the bb* rest frame. The factors 2Eb and 2Ei appear because 

the free-particle states are normahzed to 2E particles per unit volume. The factors come 

from the color traces, and the additional factor 2 in the bb case comes from the spin trace. 

From this matching process, we deduce that 

e;io _ _ „ — 

\M\^ = 2,oz.^2.oz^^2 ^'"'"'^^^'^^""^'^^ {bb\0,eS,)\bb){0\OrCSom. (11) 
Z4:6mg[Zhb) [2Ei)'' 

To compute the decay rate of the process T S'y, we replace the bb state by the T 
state and replace the bb* state by the 5* state in the squared matrix element, multiply by the 
two-body phase $2 = (l/87r)[l — (M|/M|)], and multiply by 2Mg, where Mg is the mass 
of the bb* bound state, and Mr is the T mass. ^ The result is 

^2^2^,^,2 /T-I^O ('3c MT-X /nl/nS'/l ( 



r(T^S,) ,256.-e,W^,„., „^., (T|O.C'S.)|T)(0|OfC5o)|0) 



(2M 




243 ' " ElEfrnj 

If one considers specific polarizations of the T, then the decay rate is no longer indepen- 
dent of the angle 6 between the photon and the axis that defines the direction of longitudinal 
polarization. (In T production in e~^e~ annihilation, for example, the polarization is trans- 
verse to the beam direction.) In the case of equal population of the two transverse polariza- 
tion states, dr/d{cos6) = (3/8)(l + cos^6')r, while in the case of longitudinal polarization, 
dT/d{cose) = (3/4) (1 - cos^ ^)r, where T is found in Eq. 

Using the nonrelativistic approximations Ef, ^ and E^ ^ in Eq. ()12|). we obtain 
the branching fraction 

Br(T Si) = j^.y^J^I^. X Br(T ^ ^.V^jEx,, 

81a ^ ^ mlmlm?- ^ * 

( M?\ 

X f l-^jBr(T^^V)Exp. (13) 



In deriving Eq. (fT^ we use 



r(T ^ ^v) = : ' "'.vv ' ' (14) 



M4 



^ The factor 2Mg appears for the following reason. The S state is normalized to one particle per unit volume 
in the S rest frame. In order to preserve that normalization in the T rest frame, one must multiply by 
the Lorentz-contraction factor for the volume, namely, 2M^/(2i?^), where is the energy of the S in 
the T rest frame. The factor l/(2i?^) is absorbed into the conventional definition of the phase space. 
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FIG. 2: Branching fraction for the process T S'j as a function of M^. The shaded area is 
excluded at the 90% confidence level by the T — > search of the CUSB Collaboration 

and take the vacuum-saturation approximation 

{T\0,i'S,)\T) ^ (T|V^Vx|0)(0|xVV|T), 



(15) 



which neglects terms of relative order v^. Here, and throughout this paper, v is the heavy- 
quark or heavy-squark velocity in the onium rest frame. 

In the vacuum-saturation approximation, the color-singlet production matrix element 
(0|Of (^5*0)10) may be replaced by the corresponding decay matrix element (S'|(9i(^5'o)|S'), 
with uncertainties of relative order v'^. Furthermore, the sbottomonium decay matrix ele- 
ment (S'|(9i(^S'o)|5') is related to the heavy-quarkonium (HQ) matrix element of the same 
mass by 



{s\d,cso)\s) ^ {mm's,)\m)/2, 



(16) 



where we neglect spin-dependent contributions to the HQ matrix element of relative order v"^. 
The expressions on the left and right sides of Eq. p6|) are proportional to the squares of the 
sbottomonium and quarkonium wave functions at the origin, respectively. The value of the 
HQ matrix element is known at the T mass. We estimate its value at smaller quarkonium 
masses by assuming that it scales as nig . This scaling behavior is approximately that 
which one obtains from the Martin potential \x^, which describes the J/ip and T systems 
reasonably well. Then we have 



{s\d,CSo)\s) 



\ 3/2 

nih) 
6 



(T|Oi(35i)|T)/2. 



(17) 



We use a recent lattice measurement of the bottomonium matrix element: (T|(9i(^5'i)|T) = 
4.10 ±0.42 GeV^ (Ref. Q). 

In Fig. 121 we plot the branching fraction (|13|). for several values of function of 

Mg 2mi, using the estimate of the sbottomonium matrix element in Eq. (jl7j) . Here, and 
in all further numerical estimates, we take sin^ 6^ = 1/6, rrif, = 4.73 ± 0.20 GeV, a = 1/137, 



as = 0.2 ± 0.02, and Br(T ^ At+;U-)Exp = 2.48 ± 0.06% (Ref.|l7|). This value for 
corresponds approximately to the renormalization scale rrib, which is an upper bound on the 
momentum transfer in the radiative decay process. 

In order to compare our result for the branching fraction with the experimental resolution, 
it is necessary to know the width of the S'-wave sbottomonium state. The total width into 
light hadrons is given, in leading order in a^, by the width into two gluons. This quantity 
is computed, in leading order in the nonrelativistic expansion, by Nappi jl8j : 

^ ^ 99) = ^jRm\ (18) 

where -R(O) is the S radial wave function at the origin. Using Eq. (fT7|) and taking Mg ^ 2m^, 
we have 

Att ~ 

\R{0)f = -{S\d^CSo)\S). (19) 



^ss than 10 MeV in the range 
This width is less than the 



The width of the sbottomonium state into light hadrons is 
of parameters proposed in the light-bottom-squark scenario 
energy resolution in the CUSB search for monochromatic photon signals [11]. Therefore, we 
can compare our estimate of the branching fraction Br(T —>■ S'j) directly with the CUSB 
90% confidence level for the exclusion of T — X7, which is plotted, along with our estimate, 
in Fig. El We see that, if a bound state is formed, then a part of the range of mass parameters 
proposed in the light-bottom-squark scenario is disfavored. 

There are several uncertainties in our calculation. The uncertainty in Br(T — > /i'^/i") 
is 2.4%. The uncertainty in is 4.2%. The uncertainties in the value of as and in the 
lattice computation of the bottomonium matrix element are each about 10%. There is 
also an uncertainty from the extrapolation of the HQ matrix element from Mx to M^. 
We estimate it by checking the accuracy of the extrapolation against the phenomenological 
values of the wave functions at the origin for the T and the J/ip, as determined from 
the data for the decay rates into lepton pairs combined with the next-to-leading-order QCD 
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FIG. 3: The regions of the mr-rrig parameter space that are excluded at the 90% confidence level 



The solid curve represents 



by the T X7 search of the CUSB Collaboration (shaded region) 
the central value of the theoretical calculation, and the dashed curves show the uncertainties on 
the theoretical values, as described in the text. The strip shows the region 2 < < 5.5 GeV, 
12 < rrig < 16 GeV proposed in the light-bottom-squark scenario j^. 

expressions for those decay rates. We conclude that the extrapolation error is approximately 
31% X (Mt — Mg) I (Mx — Mji^). There are uncalculated relativistic corrections of the order 
of the square of the heavy-quark or squark velocity in the onium rest frame. We estimate 
these to yield an uncertainty in the decay rate of 20% x (Mt — M^)/(Mt — Mj/^) + 10%. 
There are uncertainties from uncalculated corrections of higher order in and from the 
imprecision in the value of the renormalization scale. We estimate these by varying the scale 
of Us from mi,/2 to 2mb- This procedure yields uncertainties in of +33% and —18%. As 
Mg approaches Mx, the momentum transfer in the radiative decay becomes considerably 
less than m^, and, in this region, our choice of scale probably results in an underestimate of 
the decay rate. In this same region, we expect to find violations of the NRQCD factorization, 
which holds only for Mx — M^ ^ Aqcd- 

In Fig. El we show the region excluded by the CUSB data. We calculate the uncertainty 
band on the boundary of the excluded region by adding, in quadrature, the theoretical 
uncertainties mentioned above. We also plot the values of and nii that are suggested in 
the light-bottom-squark scenario j^. At rrig = 12 GeV, provided that the bottom squark 
lifetime is great enough to permit formation of the bb* bound state, the mass range < 
3.5toi GeV is excluded at the 90% confidence level by the CUSB data. At rrig = 16 GeV, 
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the central and upper values of the excluded range are = 3.0 and 3.6 GeV, but the 
theoretical uncertainties do not permit us to specify a lower limit at the 90% confidence 
level. One can probe the region of higher bottom-squark masses by increasing the statistics 
of the photon sample and by examining decays from bottomonium states of higher mass, 
such as the T{3S) and the T(4S'). ^ 

Bottom squarks with mass 2 GeV < < 5.5 GeV along with gluinos with mass 
12 GeV < nig < 16 GeV are proposed in Ref. ||3| to explain the larger-than-predicted 
rate for bottom quark (bb) production at the Fermilab collider. In this Letter, we show 
that CUSB data on radiative T decays already provide an important additional constraint 
on the mass ranges. The high-statistics 2002 CLEO data should either permit discovery of 
squarkonium bound states, confirm the exclusion region of the earlier CUSB data, or further 
narrow the allowed range of supersymmetry parameter space. 
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D. Son, and M. Tuts. Work in the High Energy Physics Division at Argonne National 
Laboratory is supported by the U. S. Department of Energy, Division of High Energy Physics, 
under Contract No. W-31-109-ENG-38. 
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